该内容已被发布者删除 该内容被自由微信恢复
文章于 2017年11月5日 被检测为删除。
查看原文
被用户删除
其他

同一套RNA-seq,为什么公司做的跟师兄跑的结果不一样? | TPM、read counts、RPKM/FPKM你选对了吗?

2017-09-18 小哈 嘉因生物

RNA测序分析,哪家公司适合我一文中,聊过类似的问题。


这次小哈搬来了statQuest的3个最新视频,https://statquest.org/video-index/,点击左下角“阅读原文”直达statQuest视频目录。帮你理清RNA-seq数据预处理方法,哪个更适合你的问题。


先说结论:

  1. 学术界已经不再推荐RPKM、FPKM;

  2. 比较基因的表达丰度,例如哪个基因在哪个组织里高表达,用TPM做均一化处理;

  3. 不同组间比较,找差异基因,先得到read counts,然后用DESeq2或edgeR,做均一化和差异基因筛选;如果对比某个基因的KO组和对照,推荐DESeq2。


如果找公司做RNA-seq数据处理,计算表达量时,记得要read counts。




RPKM/FPKM、TPM


https://v.qq.com/txp/iframe/player.html?vid=q05510q9xy4&width=500&height=375&auto=0

RNA-Seq分析|RPKM, FPKM, TPM, 傻傻分不清楚?里,根据上面的视频用汉语描述了这三个值的区别。


下面是视频截图

不再用RPKM/FPKM,现在推荐用TPM


一表看懂TPM更适合比较同一基因在不同sample间表达丰度的差异


DESeq2或edgeR


DESeq2和edgeR不用RPKM/FPKM或TPM做均一化,而是直接用原始的read counts做均一化处理。



DESeq2和edgeR能很好的解决这两个问题


测序深度的差异问题,用RPKM/FPKM、TPM、DESeq2和edgeR都能处理




如果组间RNA成分差异较大,怎么办?


例如liver跟spleen比,组织特异性表达基因在里面捣乱;再例如常见的某个基因KO组与对照相比,该基因成分在KO组里缺失。






DESeq2做均一化

https://v.qq.com/txp/iframe/player.html?vid=w05515xu5kz&width=500&height=375&auto=0


用edgeR作均一化

https://v.qq.com/txp/iframe/player.html?vid=g05517d3wyj&width=500&height=375&auto=0




想用ATAC、ChIP实验研究感兴趣的基因?想用已发表的ChIP-seq、ATAC-seq、eCLIP-seq、ChIA-PET、Hi-C、DNA甲基化测序、RNA-seq数据寻找线索?找嘉因生物吧!从实验、测序,到多种数据整合分析,为您一站式解决。(点击文中蓝字了解详情)




嘉因生物公众号定位:客户共性问题解答,生信学习资源导航,高通量实验导购 | 为您提供高通量实验-测序-分析-验证一站式解决方案


电话:021-61539657

Email:marketing@rainbow-genome.com
地址:上海市杨浦区赤峰路65号同济科技园1号楼611室

您可能也对以下帖子感兴趣

文章有问题?点此查看未经处理的缓存